/* TWO IMPORTANT NOTES
1) User must supply the current directory location ~10 times throughout the file. This is necessary because the clear all statements -- which purge Stata's RAM usage -- allow this file to be run from start to finish without requiring 36+GB of RAM. 
2) Synth results datasets are named in the "keep()" statements at the end of each Synth line. These contain the composition of the synthetic control, as well as the treated and synthetic time series observations. 
*/

cd ""
use replication_data, clear
gen popn=pop/100000
gen mvfatspc=farsfats/popn
sort statenum year
by statenum: gen mvfatspclag=mvfatspc[_n-1]
gen mvfatspcdif=mvfatspc-mvfatspclag
save replication_data_modelready, replace


*Synthetic California- 10 year post-intervention results horizon (Important boundary decision, the further the post-intervention results horizon, the more states that must be excluded from the donor pool) after running this model and completing the insample placebo test, I repeat this procedure for a model with a results period extending to 2015 and the donor pool restricted accordingly (i.e. many more states dropped from the matching algorithm;
clear all
cd ""
use replication_data_modelready
drop if state=="Washington"
drop if state=="Oregon"
drop if state=="Alaska"
drop if state=="Maine"
drop if state=="Hawaii"
drop if state=="Colorado"
drop if state=="Nevada"
drop if state=="Maryland"
drop if state=="Vermont"
drop if state=="Montana"
drop if state=="Rhode Island"
drop if state=="New Mexico"
drop if state=="Dist of Columbia"
synth mvfatspc mvfatspcdif(1976) mvfatspcdif(1977) mvfatspcdif(1978) mvfatspcdif(1979) mvfatspcdif(1980) mvfatspcdif(1981) mvfatspcdif(1982)  mvfatspcdif(1983) mvfatspcdif(1984) mvfatspcdif(1986) mvfatspcdif(1987) mvfatspcdif(1988) mvfatspcdif(1989) mvfatspcdif(1990) mvfatspcdif(1991) mvfatspcdif(1992) mvfatspcdif(1993) mvfatspcdif(1994) mvfatspcdif(1995) mvfatspcdif(1996) mvfatspc(1975) mvfatspc(1976) mvfatspc(1977) mvfatspc(1978) mvfatspc(1979) mvfatspc(1980) mvfatspc(1981) mvfatspc(1982)  mvfatspc(1983) mvfatspc(1984) mvfatspc(1985) mvfatspc(1986) mvfatspc(1987) mvfatspc(1988) mvfatspc(1989) mvfatspc(1990) mvfatspc(1991) mvfatspc(1992) mvfatspc(1993) mvfatspc(1994) mvfatspc(1995) mvfatspc(1996) , trunit(5) trperiod(1996) mspeperiod(1975(1)1996)  resultsperiod(1975(1)2015) unitnames(state) figure nested allopt keep(mvfatspc_2007_15_5_test)
graph save Graph mvfatspc_2007_15_CA_5_



clear all
foreach i in 1 3 4 5 7 8 10 11 13 14 15 16 17 18 19 22 23 24 25 26 28 30 31 33 34 35 36 37 39 41 42 43 44 45 47 49 50 51 {
cd ""
use replication_data_modelready
drop if state=="Washington"
drop if state=="Oregon"
drop if state=="Alaska"
drop if state=="Maine"
drop if state=="Hawaii"
drop if state=="Colorado"
drop if state=="Nevada"
drop if state=="Maryland"
drop if state=="Vermont"
drop if state=="Montana"
drop if state=="Rhode Island"
drop if state=="New Mexico"
drop if state=="Dist of Columbia"
capture noisily synth mvfatspc mvfatspcdif(1976) mvfatspcdif(1977) mvfatspcdif(1978) mvfatspcdif(1979) mvfatspcdif(1980) mvfatspcdif(1981) mvfatspcdif(1982)  mvfatspcdif(1983) mvfatspcdif(1984) mvfatspcdif(1986) mvfatspcdif(1987) mvfatspcdif(1988) mvfatspcdif(1989) mvfatspcdif(1990) mvfatspcdif(1991) mvfatspcdif(1992) mvfatspcdif(1993) mvfatspcdif(1994) mvfatspcdif(1995) mvfatspcdif(1996) mvfatspcdif(1997) mvfatspcdif(1998) mvfatspcdif(1999) mvfatspcdif(2000) mvfatspcdif(2001) mvfatspcdif(2002) mvfatspcdif(2003) mvfatspcdif(2004) mvfatspcdif(2005) mvfatspcdif(2006) mvfatspcdif(2007) mvfatspcdif(2008) mvfatspcdif(2009) mvfatspcdif(2010) mvfatspcdif(2011) mvfatspcdif(2012) mvfatspcdif(2013) mvfatspc(1975) mvfatspc(1976) mvfatspc(1977) mvfatspc(1978) mvfatspc(1979) mvfatspc(1980) mvfatspc(1981) mvfatspc(1982)  mvfatspc(1983) mvfatspc(1984) mvfatspc(1985) mvfatspc(1986) mvfatspc(1987) mvfatspc(1988) mvfatspc(1989) mvfatspc(1990) mvfatspc(1991) mvfatspc(1992) mvfatspc(1993) mvfatspc(1994) mvfatspc(1995) mvfatspc(1996) mvfatspc(1997) mvfatspc(1998) mvfatspc(1999) mvfatspc(2000) mvfatspc(2001) mvfatspc(2002) mvfatspc(2003) mvfatspc(2004) mvfatspc(2005) mvfatspc(2006) mvfatspc(2007) mvfatspc(2008) mvfatspc(2009) mvfatspc(2010) mvfatspc(2011) mvfatspc(2012) mvfatspc(2013) , trunit(`i') trperiod(1996) mspeperiod(1975(1)1996)  resultsperiod(1975(1)2015) unitnames(state)  nested allopt keep(mvfatspc_insample_`i') 
clear all
}  




foreach i in 1 4 5 7 10 11 13 15 16 17 18 19 24 25 26 28 30 34 35 36 37 39 41 42 43 44 45 47 49 50 /*51*/ {
capture noisily use mvfatspc_insample_`i', clear
sort _time
gen Ysynthetic`i'=_Y_synthetic
gen Ytreated`i'=_Y_treated
gen wweight`i'=_W_Weight
gen conum`i'=_Co_Number
keep _time Ytreated`i' Ysynthetic`i' /*wweight`i' conum`i'*/
drop if _time==.
save mvfatspc_insample_`i'_pull, replace 
}

foreach i in 1 4 5 7 10 11 13 15 16 17 18 19 24 25 26 28 30 34 35 36 37 39 41 42 43 44 45 47 49 50 /*51*/ {
merge 1:1 _time using mvfatspc_insample_`i'_pull, keepusing(Ytreated`i' Ysynthetic`i' /*wweight`i' conum`i'*/) nogenerate
}


foreach i in 1 4 5 7 10 11 13 15 16 17 18 19 24 25 26 28 30 34 35 36 37 39 41 42 43 44 45 47 49 50 /*51*/ {
gen diff`i'=Ytreated`i'-Ysynthetic`i'
gen diffsqrd`i'=diff`i'^2
egen diffsqrdmeanpre`i'=mean(diffsqrd`i') if _time<1996
egen diffsqrdmeanpost`i'=mean(diffsqrd`i') if _time>=1996
gen RMSPE_pre`i'= sqrt(diffsqrdmeanpre`i')
replace RMSPE_pre`i'= RMSPE_pre`i'[_n-1] if mi(RMSPE_pre`i') & _n>1
gen RMSPE_post`i'= sqrt(diffsqrdmeanpost`i')
gen ratio1_`i'= RMSPE_post`i'/RMSPE_pre`i'
gen ratio2_`i'= RMSPE_pre`i'/RMSPE_post`i'
}

save merged_insample_1996_10yr, replace

*CA 2015 horizon;
clear all
cd ""
use replication_data_modelready
drop if state=="Washington"
drop if state=="Oregon"
drop if state=="Alaska"
drop if state=="Maine"
drop if state=="Hawaii"
drop if state=="Colorado"
drop if state=="Nevada"
drop if state=="Maryland"
drop if state=="Vermont"
drop if state=="Montana"
drop if state=="Rhode Island"
drop if state=="New Mexico"
drop if state=="Michigan"
drop if state=="New Jersey"
drop if state=="Arizona"
drop if state=="Dist of Columbia"
drop if state=="Delaware"
drop if state=="Connecticut"
drop if state=="New Hampshire"
drop if state=="Massachusetts"
drop if state=="Illinois"
drop if state=="Minnesota"
drop if state=="New York"
synth mvfatspc mvfatspcdif(1976) mvfatspcdif(1977) mvfatspcdif(1978) mvfatspcdif(1979) mvfatspcdif(1980) mvfatspcdif(1981) mvfatspcdif(1982)  mvfatspcdif(1983) mvfatspcdif(1984) mvfatspcdif(1986) mvfatspcdif(1987) mvfatspcdif(1988) mvfatspcdif(1989) mvfatspcdif(1990) mvfatspcdif(1991) mvfatspcdif(1992) mvfatspcdif(1993) mvfatspcdif(1994) mvfatspcdif(1995) mvfatspcdif(1996) mvfatspc(1975) mvfatspc(1976) mvfatspc(1977) mvfatspc(1978) mvfatspc(1979) mvfatspc(1980) mvfatspc(1981) mvfatspc(1982)  mvfatspc(1983) mvfatspc(1984) mvfatspc(1985) mvfatspc(1986) mvfatspc(1987) mvfatspc(1988) mvfatspc(1989) mvfatspc(1990) mvfatspc(1991) mvfatspc(1992) mvfatspc(1993) mvfatspc(1994) mvfatspc(1995) mvfatspc(1996) , trunit(5) trperiod(1996) mspeperiod(1975(1)1996)  resultsperiod(1975(1)2015) unitnames(state) figure nested allopt keep(mvfatspc2015_5_dup)
graph save Graph mvfatspc_CA2015_5



clear 
foreach i in 1 4 5 10 11 13 15 16 17 18 19 25 26 28 34 35 36 37 39 41 42 43 44 45 47 49 50 51 {
cd ""
use replication_data_modelready
drop if state=="Washington"
drop if state=="Oregon"
drop if state=="Alaska"
drop if state=="Maine"
drop if state=="Hawaii"
drop if state=="Colorado"
drop if state=="Nevada"
drop if state=="Maryland"
drop if state=="Vermont"
drop if state=="Montana"
drop if state=="Rhode Island"
drop if state=="New Mexico"
drop if state=="Michigan"
drop if state=="New Jersey"
drop if state=="Arizona"
drop if state=="Dist of Columbia"
drop if state=="Delaware"
drop if state=="Connecticut"
drop if state=="New Hampshire"
drop if state=="Massachusetts"
drop if state=="Illinois"
drop if state=="Minnesota"
drop if state=="New York"
capture noisily synth mvfatspc mvfatspcdif(1976) mvfatspcdif(1977) mvfatspcdif(1978) mvfatspcdif(1979) mvfatspcdif(1980) mvfatspcdif(1981) mvfatspcdif(1982)  mvfatspcdif(1983) mvfatspcdif(1984) mvfatspcdif(1986) mvfatspcdif(1987) mvfatspcdif(1988) mvfatspcdif(1989) mvfatspcdif(1990) mvfatspcdif(1991) mvfatspcdif(1992) mvfatspcdif(1993) mvfatspcdif(1994) mvfatspcdif(1995) mvfatspcdif(1996)  mvfatspc(1975) mvfatspc(1976) mvfatspc(1977) mvfatspc(1978) mvfatspc(1979) mvfatspc(1980) mvfatspc(1981) mvfatspc(1982)  mvfatspc(1983) mvfatspc(1984) mvfatspc(1985) mvfatspc(1986) mvfatspc(1987) mvfatspc(1988) mvfatspc(1989) mvfatspc(1990) mvfatspc(1991) mvfatspc(1992) mvfatspc(1993) mvfatspc(1994) mvfatspc(1995) mvfatspc(1996) mvfatspc(1997) mvfatspc(1998) mvfatspc(1999) mvfatspc(2000) mvfatspc(2001) mvfatspc(2002) mvfatspc(2003) mvfatspc(2004) mvfatspc(2005) mvfatspc(2006)  , trunit(`i') trperiod(1996) mspeperiod(1975(1)1996)  resultsperiod(1975(1)2015) unitnames(state)  nested allopt keep(mvfatspc_insample_2015_`i') 
clear all
}  




foreach i in 1 4 5 10 11 13 15 16 17 18 19 25 26 28 34 35 36 37 /*39*/ 41 42 43 44 45 47 49 50 /*51*/ {
capture noisily use mvfatspc_insample_2015_`i', clear
sort _time
gen Ysynthetic`i'=_Y_synthetic
gen Ytreated`i'=_Y_treated
gen wweight`i'=_W_Weight
gen conum`i'=_Co_Number
keep _time Ytreated`i' Ysynthetic`i' /*wweight`i' conum`i'*/
drop if _time==.
save mvfatspc_insample_2015_`i'_pull, replace 
}

foreach i in 1 4 5 10 11 13 15 16 17 18 19 25 26 28 34 35 36 37 /*39*/ 41 42 43 44 45 47 49 50 /*51*/ {
merge 1:1 _time using mvfatspc_insample_2015_`i'_pull, keepusing(Ytreated`i' Ysynthetic`i' /*wweight`i' conum`i'*/) nogenerate
}


foreach i in 1 4 5 10 11 13 15 16 17 18 19 25 26 28 34 35 36 37 /*39*/ 41 42 43 44 45 47 49 50 /*51*/ {
gen diff`i'=Ytreated`i'-Ysynthetic`i'
gen diffsqrd`i'=diff`i'^2
egen diffsqrdmeanpre`i'=mean(diffsqrd`i') if _time<1996
egen diffsqrdmeanpost`i'=mean(diffsqrd`i') if _time>=1996
gen RMSPE_pre`i'= sqrt(diffsqrdmeanpre`i')
replace RMSPE_pre`i'= RMSPE_pre`i'[_n-1] if mi(RMSPE_pre`i') & _n>1
gen RMSPE_post`i'= sqrt(diffsqrdmeanpost`i')
gen ratio1_`i'= RMSPE_post`i'/RMSPE_pre`i'
gen ratio2_`i'= RMSPE_pre`i'/RMSPE_post`i'
}

save merged_insample_1996_2015horizon, replace

*LEAVE ONE OUT ITERATIVE POST-ESTIMATION TEST(5 models)
*drop ohio
clear all
cd ""
use replication_data_modelreadydrop if state=="Washington"
drop if state=="Oregon"
drop if state=="Alaska"
drop if state=="Maine"
drop if state=="Hawaii"
drop if state=="Colorado"
drop if state=="Nevada"
drop if state=="Maryland"
drop if state=="Vermont"
drop if state=="Montana"
drop if state=="Rhode Island"
drop if state=="New Mexico"
drop if state=="Michigan"
drop if state=="New Jersey"
drop if state=="Arizona"
drop if state=="Dist of Columbia"
drop if state=="Delaware"
drop if state=="Connecticut"
drop if state=="New Hampshire"
drop if state=="Massachusetts"
drop if state=="Illinois"
drop if state=="Minnesota"
drop if state=="New York"
drop if state=="Ohio"

synth mvfatspc mvfatspcdif(1976) mvfatspcdif(1977) mvfatspcdif(1978) mvfatspcdif(1979) mvfatspcdif(1980) mvfatspcdif(1981) mvfatspcdif(1982)  mvfatspcdif(1983) mvfatspcdif(1984) mvfatspcdif(1986) mvfatspcdif(1987) mvfatspcdif(1988) mvfatspcdif(1989) mvfatspcdif(1990) mvfatspcdif(1991) mvfatspcdif(1992) mvfatspcdif(1993) mvfatspcdif(1994) mvfatspcdif(1995) mvfatspcdif(1996) mvfatspc(1975) mvfatspc(1976) mvfatspc(1977) mvfatspc(1978) mvfatspc(1979) mvfatspc(1980) mvfatspc(1981) mvfatspc(1982)  mvfatspc(1983) mvfatspc(1984) mvfatspc(1985) mvfatspc(1986) mvfatspc(1987) mvfatspc(1988) mvfatspc(1989) mvfatspc(1990) mvfatspc(1991) mvfatspc(1992) mvfatspc(1993) mvfatspc(1994) mvfatspc(1995) mvfatspc(1996) , trunit(5) trperiod(1996) mspeperiod(1975(1)1996)  resultsperiod(1975(1)2015) unitnames(state) figure nested allopt keep(mvfatspc2015_5_drop36)
graph save Graph mvfatspc_CA2015_5_drop36


*drop ohio virginia
clear all
cd ""
use replication_data_modelready
drop if state=="Washington"
drop if state=="Oregon"
drop if state=="Alaska"
drop if state=="Maine"
drop if state=="Hawaii"
drop if state=="Colorado"
drop if state=="Nevada"
drop if state=="Maryland"
drop if state=="Vermont"
drop if state=="Montana"
drop if state=="Rhode Island"
drop if state=="New Mexico"
drop if state=="Michigan"
drop if state=="New Jersey"
drop if state=="Arizona"
drop if state=="Dist of Columbia"
drop if state=="Delaware"
drop if state=="Connecticut"
drop if state=="New Hampshire"
drop if state=="Massachusetts"
drop if state=="Illinois"
drop if state=="Minnesota"
drop if state=="New York"
drop if state=="Ohio"
drop if state=="Virginia"

synth mvfatspc mvfatspcdif(1976) mvfatspcdif(1977) mvfatspcdif(1978) mvfatspcdif(1979) mvfatspcdif(1980) mvfatspcdif(1981) mvfatspcdif(1982)  mvfatspcdif(1983) mvfatspcdif(1984) mvfatspcdif(1986) mvfatspcdif(1987) mvfatspcdif(1988) mvfatspcdif(1989) mvfatspcdif(1990) mvfatspcdif(1991) mvfatspcdif(1992) mvfatspcdif(1993) mvfatspcdif(1994) mvfatspcdif(1995) mvfatspcdif(1996) mvfatspc(1975) mvfatspc(1976) mvfatspc(1977) mvfatspc(1978) mvfatspc(1979) mvfatspc(1980) mvfatspc(1981) mvfatspc(1982)  mvfatspc(1983) mvfatspc(1984) mvfatspc(1985) mvfatspc(1986) mvfatspc(1987) mvfatspc(1988) mvfatspc(1989) mvfatspc(1990) mvfatspc(1991) mvfatspc(1992) mvfatspc(1993) mvfatspc(1994) mvfatspc(1995) mvfatspc(1996) , trunit(5) trperiod(1996) mspeperiod(1975(1)1996)  resultsperiod(1975(1)2015) unitnames(state) figure nested allopt keep(mvfatspc2015_5_drop36_47)
graph save Graph mvfatspc_CA2015_5_drop36_47



*drop ohio virginia texas
clear all
cd ""
use replication_data_modelready
drop if state=="Washington"
drop if state=="Oregon"
drop if state=="Alaska"
drop if state=="Maine"
drop if state=="Hawaii"
drop if state=="Colorado"
drop if state=="Nevada"
drop if state=="Maryland"
drop if state=="Vermont"
drop if state=="Montana"
drop if state=="Rhode Island"
drop if state=="New Mexico"
drop if state=="Michigan"
drop if state=="New Jersey"
drop if state=="Arizona"
drop if state=="Dist of Columbia"
drop if state=="Delaware"
drop if state=="Connecticut"
drop if state=="New Hampshire"
drop if state=="Massachusetts"
drop if state=="Illinois"
drop if state=="Minnesota"
drop if state=="New York"
drop if state=="Ohio"
drop if state=="Virginia"
drop if state=="Texas"

synth mvfatspc mvfatspcdif(1976) mvfatspcdif(1977) mvfatspcdif(1978) mvfatspcdif(1979) mvfatspcdif(1980) mvfatspcdif(1981) mvfatspcdif(1982)  mvfatspcdif(1983) mvfatspcdif(1984) mvfatspcdif(1986) mvfatspcdif(1987) mvfatspcdif(1988) mvfatspcdif(1989) mvfatspcdif(1990) mvfatspcdif(1991) mvfatspcdif(1992) mvfatspcdif(1993) mvfatspcdif(1994) mvfatspcdif(1995) mvfatspcdif(1996) mvfatspc(1975) mvfatspc(1976) mvfatspc(1977) mvfatspc(1978) mvfatspc(1979) mvfatspc(1980) mvfatspc(1981) mvfatspc(1982)  mvfatspc(1983) mvfatspc(1984) mvfatspc(1985) mvfatspc(1986) mvfatspc(1987) mvfatspc(1988) mvfatspc(1989) mvfatspc(1990) mvfatspc(1991) mvfatspc(1992) mvfatspc(1993) mvfatspc(1994) mvfatspc(1995) mvfatspc(1996) , trunit(5) trperiod(1996) mspeperiod(1975(1)1996)  resultsperiod(1975(1)2015) unitnames(state) figure nested allopt keep(mvfatspc2015_5_drop36_47_44)
graph save Graph mvfatspc_CA2015_5_drop36_47_44





*drop ohio virginia texas florida
clear all
cd ""
use replication_data_modelready
drop if state=="Washington"
drop if state=="Oregon"
drop if state=="Alaska"
drop if state=="Maine"
drop if state=="Hawaii"
drop if state=="Colorado"
drop if state=="Nevada"
drop if state=="Maryland"
drop if state=="Vermont"
drop if state=="Montana"
drop if state=="Rhode Island"
drop if state=="New Mexico"
drop if state=="Michigan"
drop if state=="New Jersey"
drop if state=="Arizona"
drop if state=="Dist of Columbia"
drop if state=="Delaware"
drop if state=="Connecticut"
drop if state=="New Hampshire"
drop if state=="Massachusetts"
drop if state=="Illinois"
drop if state=="Minnesota"
drop if state=="New York"
drop if state=="Ohio"
drop if state=="Virginia"
drop if state=="Texas"
drop if state=="Florida"
synth mvfatspc mvfatspcdif(1976) mvfatspcdif(1977) mvfatspcdif(1978) mvfatspcdif(1979) mvfatspcdif(1980) mvfatspcdif(1981) mvfatspcdif(1982)  mvfatspcdif(1983) mvfatspcdif(1984) mvfatspcdif(1986) mvfatspcdif(1987) mvfatspcdif(1988) mvfatspcdif(1989) mvfatspcdif(1990) mvfatspcdif(1991) mvfatspcdif(1992) mvfatspcdif(1993) mvfatspcdif(1994) mvfatspcdif(1995) mvfatspcdif(1996) mvfatspc(1975) mvfatspc(1976) mvfatspc(1977) mvfatspc(1978) mvfatspc(1979) mvfatspc(1980) mvfatspc(1981) mvfatspc(1982)  mvfatspc(1983) mvfatspc(1984) mvfatspc(1985) mvfatspc(1986) mvfatspc(1987) mvfatspc(1988) mvfatspc(1989) mvfatspc(1990) mvfatspc(1991) mvfatspc(1992) mvfatspc(1993) mvfatspc(1994) mvfatspc(1995) mvfatspc(1996) , trunit(5) trperiod(1996) mspeperiod(1975(1)1996)  resultsperiod(1975(1)2015) unitnames(state) figure nested allopt keep(mvfatspc2015_5_drop36_47_44_10)
graph save Graph mvfatspc_CA2015_5_drop36_47_44_10




*drop ohio virginia texas florida pennsylvania
clear all
cd ""
use replication_data_modelready
drop if state=="Washington"
drop if state=="Oregon"
drop if state=="Alaska"
drop if state=="Maine"
drop if state=="Hawaii"
drop if state=="Colorado"
drop if state=="Nevada"
drop if state=="Maryland"
drop if state=="Vermont"
drop if state=="Montana"
drop if state=="Rhode Island"
drop if state=="New Mexico"
drop if state=="Michigan"
drop if state=="New Jersey"
drop if state=="Arizona"
drop if state=="Dist of Columbia"
drop if state=="Delaware"
drop if state=="Connecticut"
drop if state=="New Hampshire"
drop if state=="Massachusetts"
drop if state=="Illinois"
drop if state=="Minnesota"
drop if state=="New York"
drop if state=="Ohio"
drop if state=="Virginia"
drop if state=="Texas"
drop if state=="Florida"
drop if state=="Pennsylvania"
synth mvfatspc mvfatspcdif(1976) mvfatspcdif(1977) mvfatspcdif(1978) mvfatspcdif(1979) mvfatspcdif(1980) mvfatspcdif(1981) mvfatspcdif(1982)  mvfatspcdif(1983) mvfatspcdif(1984) mvfatspcdif(1986) mvfatspcdif(1987) mvfatspcdif(1988) mvfatspcdif(1989) mvfatspcdif(1990) mvfatspcdif(1991) mvfatspcdif(1992) mvfatspcdif(1993) mvfatspcdif(1994) mvfatspcdif(1995) mvfatspcdif(1996) mvfatspc(1975) mvfatspc(1976) mvfatspc(1977) mvfatspc(1978) mvfatspc(1979) mvfatspc(1980) mvfatspc(1981) mvfatspc(1982)  mvfatspc(1983) mvfatspc(1984) mvfatspc(1985) mvfatspc(1986) mvfatspc(1987) mvfatspc(1988) mvfatspc(1989) mvfatspc(1990) mvfatspc(1991) mvfatspc(1992) mvfatspc(1993) mvfatspc(1994) mvfatspc(1995) mvfatspc(1996) , trunit(5) trperiod(1996) mspeperiod(1975(1)1996)  resultsperiod(1975(1)2015) unitnames(state) figure nested allopt keep(mvfatspc2015_5_drop36_47_44_10_39)
graph save Graph mvfatspc_CA2015_5_drop36_47_44_10_39

